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IDENTIFICATION AND MAPPING OF SINGLE NUCLEOTIDE 
POLYMORPHISMS IN THE HUMAN GENOME 

BACKGROUND OF THE INVENTION 

Field of the invention 

The invention relates to the role of genes in human diseases. More 
particularly, the invention relates to compositions and methods for identifying 
genes that are involved in human disease conditions. 

Summary of the related art 

During the past two decades, remarkable developments in molecular biology 
and genetics have produced a revolutionary growth in understanding of the 
implication of genes in human disease. Genes have been shown to be directly 
causative of certain disease states. For example, it has long been known that sickle 
cell anemia is caused by a single mutation in the human beta globin gene. In many 
other cases, genes play a role together with environmental factors and /or other 
genes to either cause disease or increase susceptibility to disease. Prominent 
examples of such conditions include the role of DNA sequence variation in ApoE in 
Alzheimer's disease, CKR5 in susceptibility to infection by HIV; Factor V in risk of 
deep venous thrombosis; MTHFR in cardiovascular disease and neural tube defects; 
p53 in HPV infection; various cytochrome p450s in drug metabolism; and HLA in 
autoimmune disease. 

Surprisingly, the genetic variations that lead to gene involvement in human 
disease are relatively small. Approximately 1% of the DNA bases which comprise 
the human genome contain polymorphisms that vary at least 1% of the time in the 
human population. The genomes of all organisms, including humans, undergo 
spontaneous mutation in the course of their continuing evolution. The majority of 
such mutations create polymorphisms, thus the mutated sequence and the initial 
sequence co-exist in the species population. However, the majority of DNA base 
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differences are functionally inconsequential in thai they neither affect the amino 
acid sequence of encoded proteins nor the expression levels of the encoded proteins. 
Some polymorphisms that lie within genes or their promoters do have a 
phenotypic effect and it is this small proportion of the genome's variation that 
5 accounts for the genetic component of all difference between individuals, e.g., 

physical appearance, disease susceptibility, disease resistance, and responsiveness to 
drug treatments. 

The relation between human genetic variability and human phenotype is a 
central theme in modern human genetic studies. The human genome comprises 
10 approximately 4 billion bases of DNA. The Human Genome Project is uncovering 
more and more of the of the consensus sequence of this genome. However, there 
remains a need to identify the nature and location of genetic variations that are 
m implicated in human disease conditions. 

Sequence variation in the human genome consists primarily of single 
15 nucleotide polymorphisms ("SNPs") with the remainder of the sequence variations 
i;S being short tandem repeats (including microsatellites), long tandem repeats 

(minisatellite) and other insertions and deletions. A SNP is a position at which two 
Q alternative bases occur at appreciable frequency (i.e. >1%) in the human population. 
^ " A SNP is said to be "allelic" in that due to the existence of the polymorphism, some 
20 members of a species may have the unmutated sequence (i.e., the original "allele*') 
whereas other members may have a mutated sequence (i.e., the variant or mutant 
allele). In the simplest case, only one mutated sequence may exist, and the 
polymorphism is said to be diallelic. The occurrence of alternative mutations can 
give rise to triallelic polymorphisms, etc. SNPs are widespread throughout the 
25 genome and SNPs that alter the function of a gene may be direct contributors to 
phenotypic variation. Due to their prevalence and widespread nature, SNPs have 
potential to be important tools for locating genes that are involved in human 
disease conditions. Wang et al, Science 280: 1077-1082 (1998), discloses a pilot study 
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in which 2,227 SNPs were mapped over a 2.3 megabase region of DNA. 

To be useful for locating and identifying genetic variations linked to human 
disease, however, it is necessary to identify and map a much larger number of Si\'P 
and to do so throughout the human genome. There is therefore a need for the 
identification and mapping of a very large number of SNPs throughout the entire 
human genome. 



BRIEF SUMMARY OF THE INVENTION 



The invention provides identification and mapping of a very large number of 
SNPs throughout the entire human genome. 

In a first aspect, the in\^ention pro\'ides SNP probes which are useful in 
classifying people according to their genetic variation. The SNP probes according to 
the invention are oligonucleotides which can discriminate between alleles of a SNP 
nucleic acid in conventional allelic discrimination assays. 

In a second aspect, the invention provides methods for using a large-scale 
map of SNPs throughout the human genome to isolate and identify genes that are 
relevant to the prevention^ causation, or treatment of human disease conditions. 
Preferred embodiments of this aspect of the invention include linkage studies in 
families, linkage disequilibrium in isolated populations, association analysis of 
patients and controls and loss-of-heterozygosity studies in tumors. 



4 



BRIRF DESCRIPTION OF THE DRAWINGS 

Figure 1 depicts the number of human restriction fragments with sizes in a 
200 bp range centered on a given point for a typical six-cutter restriction enzyme. 

Figure 2 depicts for each SEQ ID NO., the polymorphism within the 
5 consensus sequence, the position of the polymorphism in the consensus sequence 
along with the identity of the polymorphism and frequency of the alleles, and the 
map location of the identified sequence. For example, for a polymorphism in 
which "a" is identified 4 times and "t" is identified 2 times within a consensus 
sequence at position 35 from the 5' end, the text identifying the sequence will read 
□10 "SEQ ID NO. ###; polymorphism=w; position=35; alleles=a(4)t(2)." In some cases, 
the polymorphism consists of a single base deletion. In this case, the deleted base is 

1 indicated as a hyphen (-). The map location of the listed sequence is described by 
^ each of the various means which were used to identify the location, including the 

following: 

^15 1) base location relative to GenBank hit is listed as "sequence=ACC/Off" 

2 where "Acc" is the accession number of the matching GenBank entry and "Off" is 
i the offset of the polymorphism from the start of the GenBank entry, for example, 
U "sequence=M39218/98112" indicates that the polymorphism is 98,112 base pairs 

offset from the start of GenBank entry M39218. 

20 2) chromosome number is listed as chromosome=N, where N is the 

chromosome number, for example "chromosome=12". 

3) cytogenetic position is listed as cytogenetic=I, where I is the cytogenetic 
position, for example "cytogenetic=lql2.3". 

4) radiation hybrid ("rh") position relative to a GenBank entry is listed as 
25 rh=Acc/Offset (P), where "Acc" is the accession number of the relative GenBank 

entry, "Offset" is the centiray distance from the relative Genbank entry, and "(P)" is 
the radiation hybrid panel used. For example "rh=M39128/21,2 (TNG)" indicates 
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that the sequence is located 2L2 centiray from GenBank entry M39128 using the 
TNG radiation hybrid panel. Multiple map coordinates may be provided for any 
SEQ ID NO. and each coordinate is separated by a space, for example 
"map location-[chromosome-12 rh=M39128/2L2(TNG) cytogenetic-12ql8-l]." 
5 When the map position is unknown, the map fields are blank. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



The invention relates to the role of genes in human diseases. More 
particularly, the invention relates to compositions and methods for identifying 
genes that are involved in human disease conditions. Any patents and publications 
5 cited herein reflect the knowledge in this field and are hereby incorporated by 
reference in entirety. Any conflict between any reference cited herein and the 
specific teachings of this specification shall be resolved in favor of the latter. 

The invention provides identification and mapping of a very large number of 
SNPs throughout the entire human genome. This contribution allows scientists to 
10 isolate and identify genes that are relevant to the prevention, causation, or 
'S treatment of human disease conditions. 

Jj{ In a first aspect, the invention provides SNP probes which are useful in 

'"^ classifying people according to their genetic variation. The SNP probes according to 
a 15 the invention are oligonucleotides which can discriminate between alleles of a SNP 
n nucleic acid in conventional allelic discrimination assays. As used herein, a "SNP 
;:2 nucleic acid" is a nucleic acid sequence which comprises a nucleotide which is 
j=i variable within an otherwise identical nucleotide sequence between individuals or 
groups of individuals, thus existing as alleles. Such SNP nucleic acids are preferably 
20 from about 15 to about 500 nucleotides in length. The SNP nucleic acids may be part 
of a chromosome, or they may be an exact copy of a part of a chromosome, e.g., by 
amplification of such a part of a chromosome through PGR or through cloning. 

The SNP probes according to the invention are oligonucleotides that are 
complementary to a SNP nucleic acid. The term "complementary" means exactly 
25 complementary throughout the length of the oligonucleotide in the Watson and 
Crick sense of the word. In certain preferred embodiments, the oligonucleotides 
according to this aspect of the invention are complementary to one allele of the SNP 
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nucleic acid, but: not to any other allele of the SNP nucleic acid. Oligonucleotides 
according to this embodiment of the in\'ention can discriminate between alleles of 
the SNP nucleic acid in various ways. For example, under stringent hybridization 
conditions, an oligonucleotide of appropriate length will hybridize to one allele of 
5 the SNP nucleic acid, but not to any other allele of the SNP nucleic acid. (See e.g., 
Saiki et al, Proc. Natl, Acad. Sci. USA 86: 6230-6234 (1989)). For this application, 
preferred oligonucleotide lengths are from about 15 nucleotides to about 25 
nucleotides. Preferred final hybridization conditions for this application are 2x PBS 
at room temperature. Preferably, the oligonucleotide is labeled, most preferably by a 
10 radiolabel, an enzymatic label, or a fluorescent label. Alternatively, an 
i;3 oligonucleotide of appropriate length can be used as a primer for PGR, wherein the 
5 3' terminal nucleotide is complementary to one allele of the SNP nucleic acid, but 
; ^ not to any other allele. In this embodiment, the presence or absence of 
Q amplification by PGR determines the haplotype of the SNP nucleic acid. 

^^'\5 To identify the SNP nucleic acids (sometimes referred to hereafter simply as 

Q "SNPs") present in the human genome, a whole genome approach was taken to 

'p: identify SNPs on a large scale. The method described in the following examples, 

2 termed the "Reduced-Representation Shotgun" or "RRS", was utilized as it allows 

i-^^^^ the random sequencing of a specific subset (e.g., 1%) of the genome from a collection 

20 of individuals. 

Our intent was to sequence each fraction of the genomic DNA to a depth of 
2.5-5x coverage. This level of coverage was determined through a calculation of 
Poisson sampling for different levels of SNP allele frequency. Briefly, the 
proportion of SNPs identified increases with the depth of coverage of the 
25 sequencing (the sequencing of a fragment from one individual provides Ix of 

coverage and the sequencing of the same fragment from each additional individual 
provides and additional Ix of coverage), and more common SNPs are more rapidly 
detected than less common SNPs. The efficiency of detection, or number of SNPs 
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detected per ndditionnl Ix depth of coverage, however, peaks at about 2.5x coverage 
and diminishes significantly when greater than 5x coverage is obtained (calculation 
not shown). 

The distribution of restriction sites tends to be uniform across the human 
5 genome (with the exception of restriction sites containing the CpG dinucleotide). 
Thus, the proportion of the genome present in any size fraction can be varied by the 
size and extent of the fraction taken. For example, in a survey of available genomic 
sequence data on chromosomes 22 and X, the frequency and distribution of 
restriction fragments was examined, see Table 1. 

,,40 Table 1. Distribution of Restriction Fragments in Genomic Sequence, 





Enzyme 


EcoRI 


EcoRV 


BamHI 


Hindlll 


Hindlll 




Chromosome 


22 


22 


22 


22 


X 


i fs 

■si- : 


Size Range (kb) 














1-2 


40.9 


13.7 


29.7 


44.6 


67.6 


Cy 


2-3 


33 


12.6 


24.8 


32.7 


46.6 




3-4 


27 


9.4 


18.5 


26.2 


34.5 




4-5 


17.3 


9.5 


15 


20.9 


23.8 




5-7 


28.3 


15 


22.1 


25.8 


29.3 




7-9 


16.2 


8.7 


15.4 


16 


15.6 




9-11 


10 


9.1 


11.9 


8.5 


8.6 



(Values are given as number of fragments per Mb, calculated from analysis if 14Mb 
15 or 22Mb of genomic sequence on chromosomes 22 or X, respectively) 



Chromosome-specific variation of restriction site distribution is illustrated by 
a comparison of the Hindlll analysis for chromosomes 22 and X. For this reason, 
RRS plasmid libraries made using different restriction enzymes are quite useful. 
The results of restriction fragment distribution shown in Table I above indicate that 
20 for the approximately 50 Mb of chromosome 22, about 850 distinct fragments will 
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theoretically be present in a 2-2.5 kb fraction of Hindlll or EcoRI fragments, and a 5x 
coverage of the sequence of both ends of these fragments requires approximately 
11,000 reads. In practice about 25% more reads were taken as each fraction contains 
some spillover of fragments from adjacent size fractions. 

5 The number of restriction sites in the entire human genome for a typical six- 

cutter restriction enzyme can be calculated and plotted as shown in Figure 1. As 
shown in Figure 1, there are roughly 33,000 fragments in the range of 400-600 bp, and 
about 22,000 fragments in the range 1.9-2.1 kb. Each 400-600bp fragment could be 
sequenced in a single sequencing reaction, and each 1.9-2. Ikb fragment could be 
10 sequenced in two sequencing reactions, one from each end. Thus it is apparent that 
3 approximately 33,000 reads of fragment in the range 400-600bp or 44,000 sequencing 
J'jl reads would each provide Ix coverage of the SNPs present in the selected fraction of 
yn the human genome. 

The oligonucleotides according to this aspect of the invention are useful for 
'■15 identifying people according to their haplotype for a panel of SNP nucleic acids. 
m This can be acheived by obtaining a nucleic acid sample from an individual and 
!:S using the oligonucleotides according to the invention to assay for which allele the 
a individual has for a particular set of SNP nucleic acids disclosed herein, as discussed 
above. If a sufficiently large number of SNP nucleic acids are assayed, a unique 
20 haplotype can be established as a reference for that individual. Subsequently, if a 
biological sample which may be from that individual needs to be identified, e.g., for 
forensic purposes, the oligonucleotides according to the invention can be used in 
identical assays on the biological sample, and the results can be compared to the 
reference haplotype to determine whether the biological sample is from the same 
25 individual. The oligonucleotides according to the invention are also useful in 

studies to determine the relevance of various genes to the prevention, causation or 
treatment of various human disease conditions, as further discussed below. 
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Thus, in a second aspect, the invention provides methods for using a large- 
scale map of SNPs throughout the human genome to isolate and identify genes that 
are relevant to the prevention, causation, or treatment of human disease 
conditions. Preferred embodiments of this aspect of the invention include linkage 
5 studies in famiUes, linkage disequilibrium in isolated population, association 
analysis of patients and controls and loss-of-heterozygosity studies in tumors. 

The SNP map and its methods of use according to this aspect of the invention 
transform the search for susceptibility genes through the use of association studies 
and through the use of linkage disequilibrium studies. Linkage disequilibrium 
10 studies are indirect studies in which an investigator seeks to identify the presence of 
'^i common ancestral chromosomes among susceptible individuals. Association 
'JS studies are direct studies in which an investigator tests whether a genetic variant 
\n increases disease risk by comparing allele frequencies in affecteds and controls. 

Association studies make possible the identification of genes with relatively 
ini5 common variants that confer a modest or small effect on disease risk, which is 
i-i precisely the type of gene expected in the most complex disorders. Association 

studies are logistically simpler to organize and are potentially more powerful than 
m family-based linkage studies, but they have previously had the practical limitation 
'l2 that one can only test a few guesses rather than being able to systematically scan the 
20 entire genome. In the method according to the invention, association studies can be 
extended to include a systematic search through the entire list of common variants 
in the human genome to reveal the identity of the gene or genes underlying any 
phenotype not due to a rare allele. The SNP map of the human genome provided 
by the invention will make it possible to test disease susceptibility against every 
25 common variant simultaneously, for example, by genotyping a well-characterized 
clinical population with a comprehensive DNA array. 

The SNP map used in this aspect of the invention can be prepared using a 
variety of methods. One traditional method of mapping the locus of a SNP is to 
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create a PGR assay to amplify the locus and then to perform genetic mapping or 
whole-genome radiation hybrid ("RH") mapping. Another method for mapping 
the locus of a SNP is '4n silico mapping" in which the SNP and its flanking 
sequence is "BLASTed" against the publicly available sequence, such as the sequence 
5 managed by NCBI or GenBank, in order to identify the genomic overlaps that will 
positionally map the SNPs. We utilized both RH mapping and in silico mapping to 
map the locus of the SNPs. 

The location of the identified SNPs was mapped by RH mapping onto the 
existing Stanford TNG panel through developing each SNP as an STS. The TNG 
10 panel was chosen for mapping as it has been shown to order new STS's with greater 
5 than 95% confidence at 100 kb resolution. The Stanford TNG panel consists of 90 
^0 independent hybrids with an average human marker retention per hybrid of 19%. 
jjn This panel was constructed with 50,000 rad of irradiation, resulting in human 
):i chromosomal fragments 300kb average size. The practical resolution of the TNG 
Lis panel is 21 kb. One can think of the TNG panel as a "clone library", representing a 
i-i 17-fold redundancy of the human genome, with a human insert size of 300 kb and 
% 333,000 detectable ends. 

n This map can be used for conventional linkage studies in families, linkage 

disequilibrium studies in isolated population, association analysis of patients and 

20 controls and loss-of-heterozygosity studies in tumors. For example, the linkage 

disequilibrium method of Hastbacka et al, Nature Genetics 2: 204-211 (1992), can be 
used, substituting SNPs according to the invention for the RFLPs used in that 
report. Briefly, linkage disequilibrium mapping is based on the observation that 
chromosomes having a gene associated with disease which are descended from a 

25 common ancestral mutation should show a distinctive haplotype in the immediate 
vicinity of the gene, reflecting the haplotype of the ancestral chromosome. For 
example, the method is particularly useful when there is a single disease-causing 
allele with a high frequency, so that the excess of an ancestral haplotype can be 
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detected easily, and when the allele was introduced into the population sufficiently 
long ago that recombination has made the region of strongest linkage relatively 
small. Population genetics are then used to determine how much recombination 
should be expected between the gene and one or more nearby SNPs of known map 
5 location, thus locating the gene with respect to the SNP map. 



The following examples are intended to further illustrate certain preferred 
embodiments of the invention, and are not intended to be limiting in nature. 

Example 1 

OlO Cloning and identification of SNP nucleic acids 

^ Genomic DNA was isolated from a plurality of unrelated human individuals 

and approximately equal amounts from each individual was pooled. The combined 

m genomic DNA was then cut to completion with one of the following restriction 
enzymes: Hindlll, EcoRI, EcoRV, and BamHL Other restriction enzymes are also 

iF; 

0L5 useful. The digested genomic DNA was then run on a preparative agarose gel along 
n with size markers. The agarose gel containing the electrophoresed DNA was cut 
S into size fractions such that a size range of about 200 base pairs was present in each 
i^- slice (e.g., 500-700 base pairs, 1000-1200 base pairs, 2200-2400 base pairs). The DNA 
was extracted from the gel. Eluted size fractionated DNA fragments were ligated 
20 into a phosphatased vector which had been cut using the same restriction enzyme as 
was used for the digestion of the genomic DNA. Plasmid libraries were prepared by 
transforming E.coli with the ligated vectors according to well known methods of 
transformation. The plasmid libraries were tested to confirm that they contained a 
high proportion of inserts in the selected size fractionation range. 

25 Random colonies of the transformed bacteria were picked for sequencing 

from one or both ends of the genomic DNA insert. Any available method of DNA 
sequencing could be utilized, and dye terminator chemistry was preferred for its 
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optimum resolution of the heterozygotes. As the genomic DNA libraries were 
made from a pool of individuals and the DNA was size fractionated prior to 
preparation of the DNA library, each fragment in the library was sampled multiple 
times, but in almost every case each sequencing read from a given fragment is 
derived from a different DNA sample thus providing a depth of coverage of the 
DNA genomic sequences which otherwise would be unattainable. 

After sequencing of the fragments, the sequences were clustered after masking 
all known repeats. The sequences can be clustered using readily available sequence 
assembly programs, e.g. Phrap. The sequences of each cluster were compared and 
inspected for base differences, and candidate SNPs were identified at positions where 
each base was represented by a Phred quality score of >20. All sequence variants 
other than SNPs, an estimated 20-25% of the total, were also noted. All SNPs, and 
other variants, which occurred in repetitive sequences were discarded and the 
remainder were entered into a candidate SNP database. 

A subset of the candidate SNPs were verified to confirm that the majority of 
the candidate SNPs identified by sequence analysis were informative. The 
verification was done using a PGR assay to amplify DNA from several individuals, 
plus a few pools of genomic DNA from distinct ethnic groups and the PGR products 
were sequenced using dye terminator chemistry for optimum detection of 
heterozygotes. The results, not shown, of the small-scale verification indicated that 
the identified SNPs were informative. 

In this manner we were able to identify the SNPs contained within the 
specific subset of DNA which was sequenced. Through reiterative use of the RRS 
method, we were able to identify the majority of the SNPs present in the human 
genome. The identified SNPs are listed in Figure 2. 
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Example 2 

Generation of SNP maps 

Each SNP was developed into an STS and mapped using the TNG panel by 
using the method of Stewart et al. (1997) Genome Research, vol. 7, pp. 422-433. 
Briefly, oligonucleotides for PGR amplification of the fragments containing the 
SNPs were chosen using PRIMER 3.0, a software package written at the Whitehead 
Genome Center. The oligonucleotide primers were chosen according to parameters 
that generate PGR products of 100-400 base pairs in length and that allow the use of a 
single set of PGR conditions for all STSs. PGR products are assayed by ethidium 
bromide staining following agarose gel electrophoresis. An STS containing an 
identified SNP is judged successful when the primers produce a distinct PGR 
product of the expected size from total human DNA, but fails to produce a distinct 
PGR product of this size from hamster genomic DNA, In addition, each successful 
STS is PGR ampUfied on a set of approximately 90 rodent-human somatic cell 
hybrids to assure that the STS maps to a unique human chromosome. Ethidium 
stained gel images were captured using a GGD camera system and captured data was 
automatically entered into our mapping database. 

The map location for each identified SNP is listed with the SNP sequence in 
Figure 2. 

Example 3 

SNP profiling to identify an individual 

Oligonucleotides that recognize one allele of a SNP nucleic acid are 
immobilized on a filter. Preferably, the oligonucleotides comprise oligonucleotides 
complementary to at least 10 different SNP nucleic acids and are present on the filter 
in a pre-arranged array. Each filter with bound oligonucleotides is placed in 4 ml 
hybridization solution containing 5x SSPE, 0.5% NaDodS04 and 400 ng of 
streptavidin-horseradish peroxidase conjugate (SeeQuence; Eastman Kodak). PGR- 
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amplified DNA made with biotinylated primers (20 microliters) from a sample of 
blood from an individual is denatured by addition of an equal volume of 400 mM 
NaOH/10 mM EDTA and added immediately to the hybridization solution, which is 
then incubated at SS'^C for 30 minutes. The filters are briefly rinsed twice in 2x SSPE, 
5 0.1% NaDodS04 at room temperature, washed once in 2x SSPE, 0.5% NaDodS04 at 
55°C and then briefly rinsed twice in 2x PBS (Ix PBS is 137 mM NaCl/27 mM 
KCl/8mM Na2HP04/1.5mM KH2PO4, pH 7.4) at room temperature. Color 
development is performed by incubating the filters in 25-50 ml red leuco dye 
(Eastman Kodak) at room temperature for 5-10 minutes. The result is 
10 photographically recorded and the pattern can subsequently be compared with 
r\ another biological sample to determine whether the individual can be excluded as 
the source of the biological sample. 

ff^ Example 4 

::J5 Analysis of clipped reads 

i^^^^ All RRS reads were clipped of sequencing vector and low quality ends, 

Q which set a usable read length for each read. The clipped reads were 

screened for repetitive sequence with RepeatMasker, using the default 
human settings. Only reads with >=80 non-repetitive bases and >=100 
20 Phred quality (Q) >=30 bases were used in this analysis. These RRS reads 

were assembled using phrap_manyreads. Contigs with 2 or more reads must 
be aligned from a common starting point, the enzyme identified in the 
Production Protocol. High quality base discrepancies, Q>=23, were 
identified as candidate SNPs, Further restrictions on the candidate SNPs 
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were that its neighbouring 5 bases all had Q >-15, and that at least 9 
of these 10 neighbouring bases agreed with the consensus. If the number 
of detected SNPs in one clique was greater than 4 or the depth of the 
assembly (not including the genomic sequence) was greater than 5, then 
5 all SNPs were discarded for that contig. 

Example 5 
PGR confirmation of polymorphism 

5 PGR primers were designed to flank each candidate SNP, and the resulting 

I JkO fragment amplified from each of the DNAs used to construct the library. 
i;n SNPs were considered validated if at least two distinct genotypes were 
:^ observed at the candidate position (or three, if a homozygous variant was 
S observed); in addition, no position could be heterozygous in all 
Q individuals, as this would indicate a repeat sequence. 
15 

Example 6 

BLAST analysis/comparison of base call and quality 

Each sequence was blasted to a library of known repeat sequences, and any 
read containing >50% of bases in repeats was removed. The remaining reads were 
20 blasted against one another, and candidate pairs identified if they shared 
>807o sequence identity over at least 270 bases. These candidate 
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pairs were aligned using a modified Smith-Waterman alignment, and 

candidate SNPs identified (see below). Two filters were used to ensure high accuracy 
of declaring a sequence match, and to avoid inclusion of low-level repeat sequences. 
First, a pair was declared only if the sequences aligned over their entire length (save 
5 50 bp allowed on either end for sequencing end-effects), and no more than 1% of 

the bases in the alignment were candidate SNPs (see below). Second, pairs were 
then arranged into higher-order connected component groups (using 

transitivity). Component groups with more than 8 reads were removed. Paired 
sequences (see above) were run through the algorithm "SNPfinder", which 

i30 compares the base-call and quality of each position. A candidate SNP was declared if 
two basecalls were present, the Phred score of each was >20, and the 10 bases flanking 

\ the SNP (5 on either side) were of Phred quality 

:;i5 Example 7 

Cloning and sequencing to confirm polymorphism 

A pool of 10 DNAs (the Pilot Panel) or 24 DNAs (the TSC Panel) was digested with a 
restriction enzyme, size fractionated on an agarose gel, and cloned into M13-based 
vectors. Sequences were obtained on ABI 377 or 3700 sequencers. 

20 Base-calling was performed with Phrap. 
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